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Infrahumps detected in Kepler light curve of VI 504 Cygni 
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ABSTRACT 

We present a power spectral density analysis of the short cadence Kepler data for 
the cataclysmic variable V1504 Cygni. We identify three distinct periods: the orbital 
period (1.668 ± 0.006 hours), the superhump period (1.732 ± 0.010 hours), and the 
infrahump period (1.633±0.005 hours). The results are consistent with those predicted 
by the period excess-deficit relation. 
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1 INTRODUCTION 

NASA's Kepler mission has been providing high time- 
resolutio n optical photom etry since its launch on March 
6, 2009 (|Korii et alj|20ich . Designed to survey a fixed 105 
square degree field of view, it monitors the brightness of 
over 100,000 stars primarily for the purpose of detecting 
and categorizing exoplanets. However, due to the extremely 
high quality of the instrument, it has proven itself a boon 
for other astronomical communities as well, including those 
interested in the study of Cataclysmic Variables (CVs). 

CVs are classified based on their outburst properties, 
with the most active category being dwarf novae (DN). 
These systems are marked by quasi-periodic outbursts of 
varying strength, typically increasing in brightness by a few 
magnitudes. Dwarf novae are further separated into sub- 
types based on the manner in which these outbursts present 
themselves. SU UMa-type DN are characterized by the pres- 
ence of occasional superoutbursts, dramatic eruptions sev- 
eral magnitudes brighter than normal outbursts. DN that 
do not exhibit such behavior are classified as U Gem-type. 

Whether a system exhibits SU UMa or U Gem charac- 
teristics is dependent on the orbital dynamics of the system, 
and can be characterized by its mass ratio (q — M2/M1) 
and its orbital period (see Eq. [3J. Each group tends to lie 
within separate distributions on either side of a critical mass 
ratio ((/crit ~ 0.33). U Gem-type DN have long periods cor- 
responding to mass ratios greater than the critical value 
(q > 9crit) whereas SU UMa-type DN have shorter peri- 
ods corresp onding to m ass ratios less than the critical value 
(q < g crit ) (lHellier! l200ll1. 

The Kepler field of view includes 10 CVs. Of the SU 
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[/Ma-type, only V344 Lyra h as been the subject of exten- 
sive study (|Wood et al.l ' l201ll ), and to date no equivalent 
treatment has been undertaken for V1504 Cygni, which is 
the object of interest in this paper. Like V344 Lyra, V1504 
Cygni is a SU UMa-type dwarf novae, and as such features 
superoutbursts in addition to normal outburst behavior. Su- 
perhumps, periodic signals which are presumed to arise from 
the p rograde prece ssion of an elliptically elongated accretion 
disk jHellierl l200ll ). have previously been i dentified in the 
V1504 light curves and studied in detail bv lCannizzo et all 
(|201ll ). 

A thorough analysis of the Kepler lightcurve of V1504 
Cygni suggests the presence of negative superhumps (hence- 
forth referred to by the shorter term infrahumps), as well. 
These are presumed to be associated with the nodal preces- 
sion of a tilted accretion disk. 

In this paper, we outline the process of identifying and 
extracting the period of these signals in Sec. [2] and then 
show that they are consistent with those predicted by the 
period excess-deficit relation in Sec. [3] We conclude with our 
main findings in Sec. U 



2 DATA ANALYSIS 

The Kepler data used in this study are the Simple Aper- 
ture Photometry (SAP) photon fluxes recorded in the Ke- 
pler short cadence (SC) of approximately 1 minute. They are 
taken from a 550 day span starting on June 20, 2009 and 
lasting to December 22, 2010. A portion of the associated 
light curve is shown in Fig. [T] 

Gaps in the data that correspond to monthly data 
downloads occur throughout the light curve, which if not 
treated properly can distort the temporal-resolution of the 
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Figure 1. A segment of the SAP light curve for V1504 Cygni 
spanning from March 7, 2010 to June 15, 2010. Featured promi- 
nently in the center is the 5th supcroutburst observed by Kepler. 
Notable features are (1) a supcroutburst region with pronounced 
superhumps preceded by a shoulder with a profile matching a 
typical outburst, (2) outbursts that tend to gradually increase in 
magnitude, leading up to a superoutburst, (3) example data gap 
that would be corrected before periodicity analysis is performed, 
and (4) a five day window between 350 and 355 days that is used 
to generate the plots in Figures \3\ and [4] as indicated in red. 

time-series or introduce anomalous structures in the associ- 
ated power specta. These structures were ultimately elimi- 
nated through a correction process in which the gaps were 
filled with adjacent data, half from the left of each gap and 
half from the right. This has the benefit of reducing drastic 
oscillations in the high-frequency components of the power 
spectrum, enabling clearer recognition of peaks that corre- 
spond to physical frequencies of interest. 

The improved clarity does not come without cost, how- 
ever. When a longer gap (typically corresponding to sched- 
uled data transfer or unexpected downtime) appears at the 
periphery of an outburst phase, part of the outburst is often 
copied into the gap, creating the appearance of a double- 
peak. This generates an artificial low-frequency contribution 
in the power spectrum for regions containing one of these 
gaps. This manifests in the power spectrum as intermittent 
vertical banding. This banding does not affect the identifi- 
cation or extraction of periodic signals in unaffected regions, 
and is therefore ignored for the purposes of this papery 

Once the SAP light curve was extracted from the Ke- 
pler data files and corrected for relevant artifacts, residuals 
were obtained by subtraction of a boxcar-smoothed version 
of the curve. This was done largely to eliminate saturation 
in the low-frequency spectrum during subsequent analyses, 
in order to emphasize the presence of higher frequency pe- 
riodic signals associated with the orbital dynamics of the 
system. 

Following IStill et al.l (|20T0t ). the power spectrum (Fig. 
[2} was obtained from the residuals by Short-time Fourier 
transform (STFT), which here aggregates the power spec- 
tra over a 120 hour window, stepped across the duration 
of the light curve in 12 hour increments. Though there are 
no permanent periodic structures, there are three distinct 
signals that appear with varying degrees of persistence and 
intensity. The periods associated with each signal were ob- 
tained by extracting peak frequency values from the power 

1 The Kepler Science Office maintains a complete list- 
ing of known data artifacts, and regularly releases 
updates on any changes that present themselves, at 
|http: / /archie, stsci.edu/kepler/data_release. html| 



Table 1. Summary of salient results for V1504 Cyg. 



Parameter 


Symbol 


Value (units) 


Orbital period 


Po 


1.668 ±0.006 (hours) 


Superhump period 


Psh 


1.732 ±0.010 (hours) 


Infrahump period 


Pih 


1.633 ±0.005 (hours) 


Superhump period excess 




(3.84 ±0.70)% 


Infrahump period deficit 


e_ 


-(2.10 ±0.47)% 


Mass ratio 


1 


0.168 ±0.044 


Secondary mass 


M 2 


0.128 ±0.050 (M ) 


Roche lobe volume radius 


R 2 /a 


0.246 ±0.022 


Apsidal precession period 


Pap 


45.1 ±8.1 (hours) 


Nodal precession period 


Pno 


77.8 ± 17.3 (hours) 



spectrum (see Fig. [3J. In addition, the extracted periods 
were confirmed by using a Levenberg-Marquardt non-linear 
least squares algorithm to fit a sinsoidal function to selected 
quiescent regions of the light curve (see Fig. 

The most dramatic feature in Fig. [5] is the presence of 
superhumps whose contributions dominate during the super- 
outburst phase. These occur when instabilities form during 
an outburst (hence the shoulder-like feature seen proceed- 
ing the superoutburst in Fig. [T]), elongating the disk, which 
then begins to precess. These signals correspond to a period 
of P s h = 1.732±0.010 hours, in excell ent agreement with the 
superh ump period value measured bv lAntonvuk fc Pavlenkq 
|2005l ). 

There is also evidence for two additional periodic signals 
that arise in the latter half of the power spectrum. We iden- 
tify the weaker of these two signals, which becomes apparent 
approximately 200 days after the start of data acquisition, as 
the orbital period (P = 1.668±0.006 hours), in good agree- 
ment with the value measured by iThorstensen fc Taylor! 
d | l998h . Once manifested, this signal persists throughout the 
remainder of the data set. 

The second, stronger signal begins to contribute after 
approximately 250 days, and then intensifies, dominating 
the upper end of the frequency spectrum. It corresponds 
to a period of Ph = 1.633 ± 0.005 hours. There is strong 
evidence that this period corresponds to infrahumps, as will 
be discussed in Sec. [3] 



3 DISCUSSION 

It is useful to discuss the nature of a dwarf nova's orbital 
dynamics in terms of the so-called superhump period excess, 
defined as 



which is a measur e of fraction al apsidal precession. 

IPattersonl (|l998T l and lMurravl l|2000T) showed that this value 
is dependent on the mass ratio of the system, which itself 
can be recast in terms of the orbital period. This not only 
provides a useful metric for quantifying whether or not an 
observed periodic signal is associated with superhumps, but 
also for obtaining the mass-ratio due to the readily acces- 
sible measurement of orbital period. Using the formulation 
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Figure 2. The power spectrum for the region of the light curve shown in Fig. [T] highest power output appears in red whereas the lowest 
appears in blue. The superoutburst that occurs near day 310 in Fig. [T] is associated with a dramatic increase in the contribution from 
confirmed superhumps. Two additional, weaker signals are indicated by the horizontal banding observed at periods of ~1.67 and ~1.63 
hours. These are confirmed to be the orbital, and infrahump periods, respectively. 
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Figure 3. The power spectrum of the residuals for a 5 day 
window spanning from 350 to 355 (BJD-2454663) days, as indi- 
cated by the shaded box in Fig. [T] The peak at approximately 
0.01 min -1 corresponds to an infrahump period of 1.633 ± 0.005 
hours. 



presented bv lPattersonl (|2005h . the superhump period excess 
is related to the mass ratio of the binary system by 



e+ = 0.18q + 0.29g 



(2) 



With the assume d mass of a white dwa rf (Mi = 0.76 ± 
0.22 M Q ) and the lSmith fc Dhillonl (| 19981 ) secondary mass- 
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Figure 4. Example of the residual flux within the 350-355 (BJD- 
2454663) day window highlighted in Fig. [T] The resultant sinu- 
soidal fit is overlaid in blue. The fit result for the infrahump pe- 
riod is 1.629 ± 0.024 hours, in good agreement with the values 
extracted from the power spectrum. 



period relation 

M 2 « (0.038 ± 0.003)P o (1 - 58±0 09) M G 



(3) 



it is straightforward to calculate the mass ratio from an 
observed orbital period. Here, the measured superhump 
period excess for V1504 results in a mass ratio of q = 
0.1677 ± 0.0086, well below the critical value that differ- 
entiates between SU UMa- and U Gem-type systems and 
is therefore required for the emergence of superhumps, con- 
firming V1504's classification. Likewise, an estimate for the 
mass of the secondary can be placed at 0.1275 ± 0.0375 M . 

The superhump period excess is also valuable for deter- 
mining whether or not an observed periodic signal is associ- 
ated with infrahumps. Defining the infrahump period deficit 
(e_ ) in a manner identical to Eq. [T] one obtains a value re- 
lated to the fractional nodal precession (negative since the 
precession is retrograde). This has been empirically shown 
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Figure 5. Theoretical dependence on orbital period for super- 
hump period e xcess (solid line) a nd infrahump period deficit (dot- 
ted line) from IPattersonl | |2005| ) . The shaded regions associated 
with each line reflect the uncertainty in the values used in Eq. [3] 
Plotted as grey triangles a re the values for superhumps using 
data from IPatterso n] |l99 8|) and infrahumps using a selection of 
data compiled bv lWood. Thomas fc Simpson] (120091 ) from numer- 
ous sources. The V1504 values for superhump excess (red square) 
and infrahump deficit (blue square) are plotted accordingly, and 
fall comfortably within the error regions. 

by IPattersonl (|l999h to be directly related to the superhump 
period excess by e+/e_ = —2. This relationship, along with 
Eq. [2] and Eq.O is used to plot the theoretical e± values for 
varying orbital periods, along with data from several other 
systems that display super- and infrahumps in Fig. [5] The 
measured value for the third periodic signal present in the 
power spectrum resides comfortably within the error bars 
associated with e_, suggesting strongly that it is caused by 
infrahumps. Though this feature is not permanent, it does 
persist for a significant du ration after manif esting, raising 
questions as to its source. iKato et al.1 (|2012l ) suggest that 
such infrahumps may be the result of weak excitations of 
tidal instabilities in the disk that fail to trigger a superout- 
burst, which may imply a common origin for superhump and 
infrahump excitations. 

These results, along with additional calculated values 
for aps i dal an d nodal precession from formulae described in 
iHellierl (|200lh and the Roche lobe volume radius from for- 
mulae described in ISmith fc Dhillonl(| 19981) , are summarized 
in Table Q] 



4 CONCLUSION 

By examining the Kepler short cadence data for V1504 
Cygni, we have identified three distinct periodic signals. Two 
of these are readily identifiable as the orbital period and the 
well-studied superhump period, each confirmed by indepen- 
dent measurements. 

The presence of infrahumps in V1504 is inferred from 



the measurement of a third periodic signal in the power spec- 
trum that corresponds to neither the established orbital pe- 
riod nor superhump period, and confirmed by comparison of 
the infrahump period deficit against empirical models. We 
have measured the infrahump period to be 1.633 ± 0.005 
hours. 
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